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Abstract. lES 1218+304 is among a group of several 
TeV blazars that exhibit surprisingly hard energy 
spectra given their redshifts. The VERITAS collab- 
oration has carried out an intensive observing cam- 
paign of lES 1218+304 in early 2009 to characterize 
the emission properties of one of these hard-spectra 
distant blazars in more detail. While hard-spectra 
TeV blazars are most promising for unveiUng puta- 
tive signatures from EBL absorption, they may also 
have impUcations on shock acceleration scenarios in 
relativistic jets. Our observations of lES 1218+304 
demonstrate significant flux variability. Furthermore, 
we present an energy spectrum extending into the 
multi-TeV regime. The implications of these data for 
7-ray emission models are discussed. 
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I. Introduction 

With over 20 confirmed TeV blazars (see 1 1 1), TeV 
7-ray astronomy has seen an explosion in the number 
of sources over recent years. Blazars are crucial to 
our understanding of the EBL since TeV photons are 
known to be attenuated by the ambient photons via the 
77 e+e" interaction |2|[3]. Several of the newly dis- 
covered blazars at relatively large redshifts are observed 
to have hard energy spectr^. Lower limits from galaxy 
counts lfT2l suggest that the source spectra of these 
distant blazars are surprisingly hard after correcting for 
EBL absorption |13|. Consequently, models of VHE 
emission from TeV blazars that involve standard shock 
acceleration and synchrotron-self-Compton (SSC) sce- 
narios are being tested if not challenged by the presence 
of such hard intrinsic spectra lfT4ll . 

In order to better understand the emission 
characteristics of these hard spectrum blazars, the 
VERITAS collaboration has carried out an intensive 
observing campaign of the TeV blazar 1ES1218+304 
in the Spring of 2009. TeV 7-ray emission from lES 
1218+304 (redshift z = 0.182), a high-frequency peaked 
BL Lac (HBL) object, was initially discovered by the 
MAGIC collaboration in 2005 Q and soon confirmed 
by VERITAS in 2006 |8|. Previous observations 
of lES 1218+304 by VERITAS H resulted in a 

'lES 1101-232 [4], lES 0229+200 |5|, lES 1011+496 (6J, lES 
1218+304 (3, d, lES 0347-121 191, RGB J0710+591 QO), QT] 
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Fig. 1. The distribution of 9^ for excess events (ON — OFF) from 
the observation of lES 1218+304. The dashed-dotted line shows the 
boundary of the signal region (6 = 0.154°). The solid curve indicates 
the expected 9^ distribution from a point source. 



Strong detection and an integral flux measurement 
of $(>200 GeV) = 12.2 x IQ-^ photons m'^s"! 
(corresponding to ~ 5% of the flux of the Crab Nebula). 



II. Observations 

The VERITAS observatory is an array of four 12 me- 
ter imaging atmospheric Cherenkov telescopes (lACT), 
located at the F. L. Whipple Observatory in southern 
Arizona 1 15 1. Each telescope is fitted with 350 hexagonal 
mirrors following the Davies-Cotton design. The focal 
plane instrument for each telescope consists of 499 
fast photomultiplier tubes and preamps spanning a 3.5° 
field of view. The stereoscopic system combined with 
a large collection area (~ 10^ m^) and a wide field of 
view allows VERITAS to be sensitive over a range of 
energies from 100 GeV to 30 TeV. At present, the energy 
resolution of single 7-rays is 15-20% and the telescope 
array has an angular resolution of ~ 0.1°. For further 
details about VERITAS, see e.g. [23. 

lES 1218+304 was observed with VERITAS for a 
total of 39.2 hours from December 2008 through May 
2009. Data taken during bad weather, with unstable 
trigger rates or other anomalies were selectively 
removed from the data set. These quality selection cuts 
yielded a final exposure time of 30 hours at a mean 
zenith angle (ZA), of ZA=12°. All observed data were 
taken in wobble mode where the source is offset by 
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TABLE I 

Summary of observations and analysis of lES 1218+304. Total live-time, range of zenith angles for observations, 

SIGNIFICANCE OF THE EXCESS EVENTS. ALSO SHOWN, ARE THE INTEGRAL FLUX ABOVE 200 GeV ASSUMING A SPECTRAL INDEX, T, OF 
r = 3.18 AND THE CORRESPONDING FLUX IN UNITS OF INTEGRATED CRAB NEBULA FLUX OVER 200 GeV. 
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Fig. 2. The observed differential energy spectrum of lES 1218-1-304 
from the 2008-09 data set (preliminary). Solid line shows the best fit 
power law to the data (filled circles) with a spectral index value of 
— 3.18±0.09i,tat . This agrees, within statistical and systematic uncer- 
tainties, to the previous VERITAS measurement of — 3.08 ± 0.34stat 
(dashed Une) (t). 



0.5° relative to the center of the field of view (FoV). 
This allows us to simultaneously sample background 
events, utilizing part of the camera that is identically 
offset from the center of the FoV but away from the 
source. 
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Fig. 3. The daily light curve in units of integral flux above 200 GeV, 
<j>{¥. > 200GeV) from lES 1218-H304 assuming a spectral shape, 
dN/dE oc E"'" with F = 3.18. The inset shows the flaring nights 
in more detail. 



source region of radius 6 ~ 0.154° centered around 
the direction of lES 1218H-304 is used to calculate 
the ON counts. The background or OFF counts are 
estimated with the Reflected — Region 1191 model 
using seven off-source regions. The set of 7/hadron 
separation cuts for this work was optimized for a source 
with a 10% Crab-like flux. The significance of any 
excess {ON — OFF) events is calculated following the 
mefliod of Eq. (17) in Li & Ma EqI. 



III. Data Analysis 

Following calibration of the camera gains with 
nightly laser runs, the data were processed with 
several independent analysis packages. During the 
analysis process, the shower images are cleaned and 
characterized with the standard Hillas parameters [17] 
by computing various moments of the light distributions. 
Subsequently, stereo images from different telescopes 
are combined to reconstruct the arrival directions 
on the sky and the impact point on the ground. 
Gamma-ray like events are identified using cuts on the 
mean-reduced-scaled length and mean- reduced-scaled 
width parameters (see i 1 8J for a description of the 
parameters). The energy of reconstructed events is 
estimated using Monte Carlo simulations. A circular 



IV. Results 

The result of VERITAS observation and analysis of 
lES 1218H-304 in 2008-09 is summarized in Table U 
along with the results from 2006-07. A significant excess 
of 23.5(7 is detected from the direction lES 1218H-304. 
Figure [T] shows the distribution of the squared angu- 
lar distance, 0^, between the reconstructed and source 
position for any excess events {ON — OFF). The 
distribution of excess events (Fig. [U clearly peaks at 
small values of 0^, consistent with the detection of a 
point-like source. 

The time-averaged differential energy spectrum ob- 
served from lES 1218H-304 during the 2008-09 observ- 
ing season is shown in Fig. |2l The energy spectrum 
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Fig. 4. The yearly light curve in units of integral flux above 200 GeV, (pCE > 200GeV) from lES 1218+304 is shown (binned by day). The 
filled circles coiTespond to VERITAS data from the 2006-07 epoch. The open circles represents the most recent data from the 2008-09 season. 



ranges from 200 GeV to 1.8 TeV. The best power law 
fit to the measured spectral points yields, 



dN E 
dE ^0.5TeV 



-3.18±0.09a 



(1) 



where the flux normalization constant at 500 GeV is 
(10.94 ± 0.68) • 10-^ TcV-im-2s-i with a of 3.4 
for 7 degrees of freedom. 

Figure [3] shows the daily integrated photon flux above 
200 GeV for the 2008-09 observing season. The pho- 
ton flux was calculated assuming a differential energy 
spectrum of dN/dE cx E"'^ with r=3.18. The average 
integral photon flux above 200 GeV is found to be 
(nA±0.89stat) X 10~*^ m"2s-i corresponding to -7% 
of the Crab Nebula's flux. Strong nightly flux variations 
are evident in the daily light curve of lES 1218-1-304 
between January 25 and February 5, 2009. During the 
highest state, the 7-ray flux reached '^20% of the Crab 
Nebula flux. The chi-squared probability of constant 
emission during the entire season (Fig. [3]- dashed-dotted 
line) is 2 • 10~^^. The inset within Fig. [3] depicts the 
flaring activity in greater detail. The highest flux seen 
from lES 1218H-304 occurred on MJD 54861 (January 
30, 2009). The high state data in January of 2009 
indicates variability (flux doubling) on a time scale of 
days. 

Figure |4] shows the daily integral photon flux for both 
the 2006-07 and 2008-09 seasons. We note that, except 
for the flaring states, the 7-ray emission is relatively 
constant over both periods of observation. 

V. Discussion & Conclusions 

The VERITAS observation of lES 1218H-304 shows 
that the observed energy spectrum is considerably hard 
given its redshift. Moreover, the light curve clearly 
shows evidence of strong flux variability in the VHE 
gamma-ray emission of lES 1218-1-304. 



Bottcher et al. 1141 recently suggested a possible 
explanation for the hardness of the VHE spectra of high 
redshift 7-ray blazars. According to this model, TeV 
blazars may contain a separate but very hard, slowly 
varying emission component beyond several hundred 
GeV. In the presence of a low magnetic field, electrons 
in the extended kiloparsec scale jets of TeV blazars 
may inverse-Compton scatter ambient cosmic microwave 
background (CMB) photons to TeV energies. Given the 
required radiative cooling time scale of the process, this 
TeV inverse-Compton emission is predicted to be very 
slowly variable in a low B-field. 

Considering the time scale of the observed flux vari- 
ability from lES 1218-1-304, we can constrain the size 
of the emission region, R, of the VHE photons to 
i? < c • At • (5 w lO^^S pc, where the S is the relativistic 
Doppler factor of the radiating region. For a resonable 
value of the Doppler factor (<90), R ^ kpc. Therefore, 
we can rule out the CMB-Inverse-Compton model of 
emission in the kpc scale extended jets as the sole 
explanation for extreme hardness in the intrinsic spectra 
of TeV blazars. 
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